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X-ray Image Processing 

Introduction 

The present invention relates to improvements in the processing of digitised x-ray 
images (particularly of the breast - termed mammograms or mammographic images), and 
more particularly to the enhancement of such images to assist clinicians in making accurate 
diagnoses based on them. 

X-ray imaging is used as a basis for many medical techniques and, in particular, 
mammography continues to be the examination of choice for early detection of breast cancer 
in post-menopausal women and is the basis for national screening programmes. 

Early detection of breast cancer greatly improves mortality rates, perhaps by as much 
as 25%. For this reason, mammographic examinations are nov^adays performed on about 25 
million women annually in the EC (of which, about 3 million are in the UK), at a cost of 
about 3 Bn US dollars per year. This huge cost and the poor accuracy in diagnosis (8-25% of 
cancers are missed and 70-80% of open surgical biopsies turn out to be benign), have led to 
increased interest in applying computer-aided techniques. Further, pressure for a reduction in 
the screening interval and for routine adoption of two-view screening (cranio-caudal and 45 
degree medio-lateral) would entail a large increase in the number of mammograms to be 
analysed in the screening programme and this makes the development of reliable and robust 
computer techniques vital. 

As an application of image processing, mammographic images pose a tough challenge 
because they have poor signal-to-noise ratio. This is largely because the images exhibit 
complex textures, and because there is a compromise between radiation dose and image 
quality. Worse, abnormalities appear as quite subtle, irregular, often non-local differences in 
intensity and the images are inevitably cluttered due to superimposition in the image of 
features separated in the breast. Further, the background varies greatly between different 
breasts, and there is relatively weak control of the imaging acquisition. 

Unfortunately, while there have been proposals for the application of image 
processing to mammography, the vast majority have been of limited scope and incorporate 
only general non-mammography-specific image processing considerations. This involves 



great dangers. Image smoothmg using such standard techniques may make lesions easier to 
locate, but can remove calcifications and spiculations which can be signs of cancer Ed^e 
sharpening using standard techniques may appear to improve an image, but malignant lesions 
typically have fuzzy edges whereas benign ones tend to have sharp edges, so this edge 
5 sharpening process can actually transform an image of a malignant lesion into one that 
appears to a radiologist to be benign. 

Calcifications present particular interest and problems. Localised cancers in the ducts 
or lobules of the breast are often are associated with secretions that thicken or become 
necrotic. These are called calcifications, or microcalcifications if they are smaller. Clusters of 
10 microcalcifications, which appear as small bright white objects in a mammogram, can be one 
of the earliest signs of breast cancer. Thus the identification of calcifications is a major goal 
of screening programmes, though benign calcifications are common (for example they often 
occur in blood vessels), and so the distinction between ductal and vascular 
microcalcifications needs to be made if the number of false positives is to be sufficiently low 
15 This is not a problem for radiologists interpreting mammograms against the background of 
their knowledge of breast anatomy. It is. however, a challenge for image analysis programs 
Further, dust and dirt entering the imaging system can create artifacts that mimic the 
appearance of microcalcifications and thus can be a cause of false positives both for 
radiologists and automated imaging systems. 
20 Further problems are caused because the imaging process itself introduces a number 

of variables which affect the image and these will be explained below. 

In the accompanying drawings Figure 1 shows a schematic representation of the 
components ofa conventional screen-film mammographic system. When a mammogram is 
perfonned. a beam of x-ray photons I from an x-ray tube 2 powered by a generator 4 is 
25 directed towards a breast 3 compressed between compression plates 17. This beam 3 is 

filtered by filter 5 to remove low energy photons and collimated by collimator 7 to the area 
of interest. The beam has a spectrxim of energies that is characteristic of the tube voltage and 
in particular, the material of the anode 9 but the spectrum is independent of the woman being 
scanned and the view taken. The intensity of the beam exiting the breast is related to the 
0 thickness and type of tissue in the breast. The x-ray photons leaving the breast nonnally have 
to pass through an anti-scatter grid 1 1 before they reach a phosphorous intensifying screen 



13. If an x-ray photon is absorbed in the screen 13, light photons are emitted by the phosphor 
and these light photons expose a film 1 5 which is processed to produce a mammogram. The 
exposure to the breast is stopped once an automatic exposure control 19, positioned under a 
section of the breast, has received a set exposure. To generate a digital image the x-ray film 
5 is typically digitised using a laser scanner system or CCD and light box (not illustrated). 

The intensity of radiation incident on the breast in such a system varies spatially for 
several reasons. The most significant is the "anode heel effect". An x-ray tube produces 
x-rays by firing an electron beam at an anode. As the electron beam penetrates the anode the 
electrons are absorbed at varying depths and the x-ray photons that are produced have to 
10 travel through varying thicknesses of anode material before leaving the anode. This leads to 
varying attenuation of the emergent x-ray beam thus giving spatial variations in the incident 
x-ray spectrum, this is termed the anode heel effect and is quite substantial. Another source of 
spatial variation is due to the diverging nature of the beam. This means that the further away 
from the source the more spread out the x-ray beam is. However, this effect is small given 
5 that the distance from the source to the breast is large relative to the breast size. In visual 
assessment of mammograms the clinician mostly considers local variations in intensity, so 
the smooth change caused by the anode heel effect is not too troublesome. However, it does 
cause a problem for automated systems. 

Two further effects of the image forming process which affect the image are scatter 
and extra-focal radiation. Considering scatter first, the x-ray radiation passes through the 
breast as shown in Figure 2 and what is known as the primary beam passes through in a 
straight line from the anode to the intensifying screen. However, some x-ray photons are 
scattered in the breast and arrive at the screen from unexpected directions. The anti-scatter 
grid, which typically consists of a series of angled lead strips 21 separated by paper and 
aligned with the primary beam, removes many but not all of these scattered photons. Thus 
some scatter will reach the screen and be recorded on the film. Scatter can be estimated using 
the techniques described in "Computing the scatter component of mammographic images", 
by Ralph Highnam, Michael Brady and Basil Shepstone published in EEEE Med. Imaging, 
1994, 13, pp 301-3 13. This allows the calculation of the primary energy imparted to the 
screen to be improved by the subtraction of the energy due to scatter. Extra-focal radiation 
refers to radiation which comes around the edge of the collimator 7 as shown by numeral lA 




in Figure 3, and can constitute up to 15% of the total, some of which is scattered and reaches 
the film/screen. Simple techniques for measuring the extra-focal radiation are known, for 
instance from the paper by Highnam, Brady and Shepstone mentioned above. 

The processes of intensification by the intensifying screen introduces blur or glare 
5 into the recorded image because the absorption of an x-ray photon at a point site 13 A in the 
screen results in the approximately isotropic emission of light as shown in Figure 4 which 
results in blurring of the image recorded on the fihn. 

Further, the relationship between the density of the image on the film and the energy 
imparted to the intensifying screen is not linear and changes with film processing conditions. 
1 0 Again, this may not affect visual assessment which is based on local variations, but would 
affect automated analysis, especially if the non linearity were not explicitly taken into 
account. 

Finally, the process of digitizing the film introduces digitizer blur into the digital 
representation of the image. 

1 5 It will be appreciated, therefore that the enhancement of x-ray mammograms and 

automated recognition and differentiation of features in mammograms is a very difficult 
problem. 

It has been proposed that the conversion of a digitised mammogram into a particular 
representation, ternied the K, representation, is capable of improving the enhancement and 
20 analysis of such mammograms. This was described in "A representation for mammographic 
image processing" by Ralph Highnam. Michael Brady and Basil Shepstone published in 
Medical Image Analysis; 1 996, vol. 1, no. 1, ppl-18 and, since the present invention is 
concerned with improvements to it, will be briefly explained below. 

The intensity of a mammogram at a given pixel (x, y) indicates the amount of 
25 attenuation (absorption and scattering) of x-rays in the pencil of breast tissue vertically above 
(x, y) on the film. 

Ideally, one might hope to be able to produce a quantitative three-dimensional 
representation of the breast with each voxel labelled with a tissue type, such as: glandular, 
fibrous, cancerous, fat, calcium. Given the x-ray attenuation within a voxel it is certainly 
30 possible to classify fat since it has relatively low linear attenuation coefficients. It is also 

possible to classify likely occurrence of calcium, which is practically radio-opaque. However. 
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the remaining breast tissues are those that comprise anatomically significant events in breast 
disease, such as cysts, malignant masses, fibroadenomas, and they are difficult to resolve 
from x-ray attenuation measurements alone. In the h,„, representation these remaining tissues 
are classified as "interesting tissue". Further, there is actually very little calcium so for 
5 practical purposes it can be ignored. 

Unfortunately, a further problem arises because of the projective nature of 
mammographic imaging: the three-dimensional information is lost. In light of this, the only 
information that is available describes the tissue within a cone of the breast, where the cone 
has as its base the area of a pixel and as its apex the x-ray source. After appropriate correction 
1 0 the x-ray beam within this cone can be considered as a pencil beam. Thus in the /?,„, 

representation (with calcium ignored) there are basically only two tissue classes of fat and 
interesting tissue to consider, and the thicknesses of the interesting tissue (/»,„, cm) and fat 
(hfa, cm) which together must necessarily add up to the total breast thickness H ( i.e. H = hi„ + 
hfc^ are used as quantitative breast measurements. 
^ ^ ^'^ practice /r,„, is computed from a mammographic image using data related to system 

calibration and image calibration. The x-ray tube output spectrum is assumed to be relatively 
stable but the anode heel effect is corrected for. The mammographic imaging process has 
several parts which might vary from day to day. In order to effect meaningful image analysis 
by computer, it is necessary to know these variations in order to make the images conform to 
20 a standard. To achieve this requires calibration data. The film-screen response, film processor 
and film digitizer are calibrated by collecting the following data; 

1 . A step wedge film. A film is produced with a stepped wedge made of lucite 
placed along the back of the film and a lucite block placed over the automatic exposure 
control. This fihn allows us to calibrate the film-screen system and fihn-processing so that 

25 energy imparted to the intensifying screen can be related to film density. 

2. A "blank"' film: A film is taken with a short time of exposure with no object 
(breast) present. The exposure has to be short so that the film does not saturate. An exposure 
of 0.04 seconds, at 100mA and 28kV for example produces a film that has film densities that 
vary between 1 .8 and 2.6 (despite looking black). This film provides information about the 

30 spatial variations of the incident radiation intensity. 
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The digitized image of step wedge film: The film density on each step of the 
wedge is measured so that once digitized, the relationship between pixel value in the digital 
image and film density in the corresponding area of the film is known. 

As well as calibrating the system components, data specific to each mammographic 
examination is needed. In particular: 
The tube voltage (V^,, kV); 
The tube current (4^, mA); 
The time of exposure (/^s); 
The breast thickness (Hem). 

Most of this information is readily available but measuring the breast thickness H is 
currently awkward since the radiographer has to measure it using a ruler, though newer 
machines are incorporating automatic measurement of breast thickness. 

Given a mammographic image, the thicknesses of interesting and fatty tissue between 
the x-ray source and each pixel can be found by considering the energy imparted to the 
intensifying screen at each pixel which is obtained from the pixel values in the image using 
the calibration data. Let E,,/x. y) be the energy imparted to the screen in the area 
corresponding to the pixel (x, y). E,„(x, y) contains both scatter and primary components. 
The primary component y) is determined by subtracting a scatter estimate from the total 
energy imparted as mentioned above. 

Now for a pixel with h,„. cm of interesting tissue and V cm of fatty tissue above the 
corresponding area of the intensifying screen, the total attenuation at any energy E is 
expected to be: " 

{E,x,y) = h im(x, y)n s^R) + hfa,{x, y) ^fa,{ E) 

= h inix,y){fi i„t(^) - + HflfariE), (1) 

25 where the substitution h^Jx. y) = H - h,Jx, y) is made. 

In this case, the energy expected to be imparted to the imensifying screen by the 
primary photons is: 

Ep{x, y) = (t>(y,ubc, X, y)Apts / N^''{E)ES{E)G{E) x e-'^'^<^)'^'-g-w,x,y)^ 



20 
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where (z) is the photon flux for an x-ray tube voltage of F^j.. this varies across the image due 
to the anode heel effect; Ap is the pixel area; is the time of exposure; No"'(E) is the 
relative number of photons at energy £; S(E) is the absorption ratio of the screen to primary 
photons of energy £; G(E) is the transmission ratio of the grid for primary photons of energy 
E; n,^(E) is the linear attenuation coefficient of the (typically) lucite (compression plate) at 
energy E\ and hp,„, is the thickness of the compression plate, all of which are known from the 
calibration or image conditions. 

Note that after substituting Equation (1) into Equation (2) the only unknown is hfjx, 
y). This can be found by equating the primary energy found in the practical case (i.e. 
measured from the image) with the theoretical value (i.e. the expected value calculated above) 
and solving the resulting nonlinear equation to determine h,Jx. y). 

This process of converting the image into the /?,„, representation can be visualised as 
converting the original image so that the fat has risen to float on top of the interesting tissue 
surface, then the fat is peeled off leaving the representation h^Jx, y). Informally, this 
15 representation can be viewed as a surface and clinically significant effects such as masses 
appear as features on this surface, eg. small hills, as seen for example in Figure 5. Note that 
this is fundamentally different from regarding the intensity image as a surface, since the A,„, 
representation is a quantitative measure of anatomical tissue which is distributed through 
vertical pencils of the breast. The importance of A,„, stems from the fact that it factors out the 
20 imaging parameters particular to the examination to yield a representation of the intrinsic 
anatomy that is ultimately what is relevant for diagnosis. 

Summary of the invention 



The present invention is concerned with improvements in the calculation of the h,„, 
representation, which not only improve the accuracy of the A,„, representation but also provide 
some very useful results along the way. 

A first aspect of the invention provides a method of correcting a digital representation 
of an x-ray image for degradation caused by a digitiser, the digital representation being a 
measurement of the image density for each of a plurality of pixels of the image, the 
measurement being obtained by illuminating the x-ray image, measuring the attenuation of 
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light by tf,e x-.y taage and calculating from the attenuated light values of U,e ima.e 



density, 

the method comprising the steps of: 



calculating fro. the pixe, values in the digital representation the value of the intensity 
of said attenuated light corresponding to each pixel; 

applying to the intensity values a modulation transfer function to correct for the 
degradation introduced by the digitiser; and 

converting the corrected intensity values back into values representative of the inaa^e 



intensity. 

10 A second aspect of the invendon provides, in a method of x-ray imaging usin. an 

.ntensrfying screen to receive x-rays and emit light U> he recorded on an x-ra'y film, a ;ethod 
of calculattng from the x-ray image density the energy which was impa„ed to the intensifying 
screen, comprising the steps of 

calibrating the x-ray fita and intensifying sc«en by measuring fl.e response of the x- 
ray film and mtensifying screen R, a plurality of differem intensities of received x-rays 
fmmg a theoretical model of the expected response to the measured response ' 
and using the fitted theoretical model to calculate the imparted energy from the x-ray 

image density. ^ 

A third aspect of U.e invention ptovides, in a meU,od of x-,.y imaging in which an 
.ntenstfymg screen is used to receive x-rays and emit light to be recorded on an x-ray fiim 
me .mage recorded on the flta is digidsed to give a pixelised digital representaUon of the ' 
m,age density on the fita. a method of enhancing the digital representation of the image to 
remove the contribution thereto of glare from the intensifying screen comprising .be steps of 

convertmg the digital representation of image density into a representation of the 
energy miparted to the screen; and 

deconvolving the representation of the energy imparted ,o the screen using a 
wetghting mask defming the point spread function for the intensifying screen remove the 
contribution thereto of glare from the intensifying screen. 

A four* aspect of the invenuon provides a method of calculating from a 
mammogram the compressed thiclcness of the imaged breas, comprising «,e step of 
dehmtting in the mammogram the region corresponding «, the part of the breast which is 
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compressed from the region corresponding to the uncompressed breast edge by detecting the 
smoothness of curves of equal intensity in the mammogram. 

The aspects above may be combined with other known steps to improve the 



production of an /7,„, representation of the breast. 

Another aspect of the invention provides a method of detecting microcalcifications in 
a breast from a mammogram of the breast, comprising the steps of processing the 
mammogram to produce an /?,„, representation according to the above methods, converting the 
value oih,„, for a candidate region mto a value representative of the volume of interesting 
tissue and thresholding the volume values to detect as microcalcifications areas of the 



^ 1 0 mammogram where the volume value exceeds a threshold. 

The invention also provides apparatus for carrying out the above methods. Further, 
the invention also provides a computer program, which can be provided on a computer- 
readable storage medium, for controlling a computer to carry-out the methods. 

15 List of drawings 

The present invention will be further described by way of example with reference to 
the accompanying drawings in which :- 



Figure 1 is a schematic representation of a conventional screen-film mammographic 



Figure 3 is a schematic view of a collimator showing the path of extra-focal radiation; 

Figure 4 is a schematic view of an intensifying screen and film; 

Figure 5 shows an h,„, surface of a mammogram; 

Figure 6 shows a film-screen characteristic curve; 

Figure 7 shows an example of fitting the characteristic curve; 

Figure 8 shows an example of the effects of compensating a mammogram for 

intensifying screen glare; 

Figure 9 shows two film density profiles with confirmed shot noise in them; 

Figure 10 shows an original energy profile and the energies after glare compensation; 

Figure 1 1 shows an example of noise detection in a phantom; 

Figure 12 shows an example of noise detection on a real mammogram; 



20 



system; 




Figure 2 is a schematic view of a typical anti-scatter grid; 
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F.gure 13 shoes a funher example of noise detection on a real .an.n.ogran,- 
Figure 14 illustrates the detection of low contrast calcifications- 
Figure 15 illustrates an image sample containing shot noise- 
Figure 16 schematically illustrates a cross-section across a compressed breast- 
Figure 17 is a schematic of a cranio-caudal mammogram; 
Figure 1 8 shows examples of the breast edge; 
Figure 19 shows examples of the breast edge with varying H- 

Figures 20 and 21 illustrate roughness measures used in L Jstimation of breast 

thickness; 

Figure 22 shows A. cooniina^ system used for e«ra-focal con,pe„sa,io„ 

Figure 24 shows «vo examples of *e extra-focal component of the radiation 

Ftgure 25 shows the modulation n^fer functions of some film digitize.. 
. Ftgure 26 illustrates the modulation t^sfer Wtions of intensifying sereins for 

dtfferent values of the attenuation coefHcien, for light in phosphor 

Ftgure 27 illustrates the modulation transfer functions of intensifying screens for 
d.,ferent values of the attenuation coemcient forx-rays in phosphor- 

thiCessT """^ ''''^ "--""^ ''"^-■« 

differenlriir^''^""-*-^'''^-^^^ 

Figure 3 1 illustrates the modulation transfer Wions of two intensifymg screens for 
different gaps between the film and screen; and greens for 

Figure 32 illustrates the geometry for calculation of the point sptead fitnctions. 

Description of preferred embodiments 

An embodiment of the invention will now be described by going through U.e steps 
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necessary to calculate the representation. These can be summarized as foUows:- 

(1) Convert pixel value P(x. y) to film density D(x. y) using the digitiser 
calibration data; 

(2) Remove digitiser blur, 

5 (3) Convert film density D(x, y) to energy imparted to intensifying screen E'-^^.fx, 

y) using film-screen calibration data; 

(4) Compensate £""Vrx, j/^ for intensifying screen glare; 

(5) Compensate E-^Jx. y) for the anode-heel effect and diverging x-ray beam; 

(6) Estimate the thickness of the compressed breast; 
^ ^ O) Estimate the scattered radiation E'"''/x, y); 

(8) Estimate the extra-focal radiation E'"''/x, y) components; 

(9) Compute the measured primary energy E^^p. 

( 1 0) Compare this measured primary energy with the theoretical primary energy 
calculated by equation (2) for different values of h^/x. y). Equating the measured and 
theoretical values using a numerical approach solves for hi„/x, y). 

Digitisation 



15 



20 



25 



Considering steps (1) and (2), mammographic images are digitised by illuminating the 
mammogram with light and measuring the amount of light transmitted through it. Various 
arrangements are used such as laser-scanning devices or devices which use a light box to 
illuminate the mammogram and a CCD camera to record the transmitted light. Digitisers are 
designed to record and output the pixel values related to the fihn (image) density rather than 
the transmitted light itself because the film density is independent of the illumination level. 

Step (1) of the method above is relatively straightforward because modem high- 
quality laser scanning devices, for instance, have a known linear relationship between the 
film density and the pixel value: 

P(^x,y) = mD(x ,y)+c, where m and c are constants. 

30 Thus D(x, y) can easily be found from the pixel value. 
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or ca„ ,e.ea.u.=d. For instance, o,odu,a.on .ansfe. <U„«ions for a 
».c™de„..„«., ceo ca„e. and a .as. scanner are disclosed in -Digi., Ma^^og^p,, 

: ~ :::r: ^ ^ " ' — 

uu.yju yjj, iyy3 and are shown in Fieure 25 VJ\fh fU^ 
functions are applied no, ,o ,he fUn, d , 

related u> the n H ' '° "W-^h 

related a the film density in a non-linear way as follows:- 

//(j:,>')x lO""''''' 

m "8^' -'-'-ed and the ntodnlation transfer functions L 

*e d,g.t.e. are applied to it to remove the digitizer hlur. Then the transmitted Hgh, is 

Film screen response 



The next step, step (3) is to calculate from the dehlutred film density D^. y, the 

The film-seen response to energy impaned to the h,.e„sifying screen is given by a ' 
W«e„s..c curve ,i.ethatinPi^.e«.Pign.e.showsafilm-screen 
^™d hy exposmg a lucite step wedge. Lncite is chosen because of its similar ahsotption and 
scattenng properttes to a hteast consisting of half fa, half interesting tissue. A small 

~ for scanered tadiation can be made, and the logarithm of the calculated relative 
energy „ ,s p,„„., „^ ^^^^ ^^^^ ^ ^ 

Wtth the mvention a setpentine characteHstic cu,ve is fitted to the lucite step wedge 
dat. and thrs fitted relationship is used to calculate the imparted enetgy fiom the fiL 
denstty. In more detail, the serpentine curve can be expressed as: 

rearranged this is: 
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y=~2 — r 



This has the properties: 



Intercept at (0, 0) 
Extremaat (±a,± ") 



Inflection at x = 0,x = ± 



-b^ b^ 
Extent: -inf < x < +inf, < V < 

2a 2a 

Symmetry (0,0) - rotational 
Asymptote y=0 

Oradien, 



Gradient at origin — = — ^ 



. ^_b^ 
dx a' 
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Thus a,b and also the origin to use need to be detennined. The origin has to be the 
place of symmetry so there are four unknown parameters: a. b. x^^ y^^^. The equation of 

the curve becomes: 

y = yongin + " ^''"Si'^) 

(jC — Xorigin)^ 

Where x = log £ (where £ is the energy imparted to the intensifying screen) zndy = D 
(where D is the film density). The origin is initially guessed to be in the centre of the 
film-density range and the film density range is measured so as to determine initial estimates 
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of a,b. An optimizing routine is used to find the optimal values of the parameters based upon 
the step wedge and blank film calibration data collected from the mammography unit. 

Inverting gives: 

...... . b^i^b'-Aa\y-y„,^Y 

— y origin 

Figure 7 shows an example of the curve fitting. In Figure 7 the horizontal axis is the 
logarithm of the energy imparted to the intensifying screen. The vertical axis plots the film 
densities. The line marked step wedge data is the most accurate and is based upon 
measurements with a piece-wise linear fit between the points. The short dashed line is the 
serpentme fit to that data showing a good fit. BCU is data collected from the Oxford Breast 
Care Unit and the longer dashed line then is the serpentine fit to the data collected at the BCU 
which gives a good fit in the low film densities, but not such a good fit at high film densities 
The fit can be improved to give the curve with long and short dashes by noting that although 
the fihn-screen curve is not symmetric in its entirety, the highest film densities are not used, 
so D„„ can be curtailed so that symmetry is obtained. 

Glare removal 



) 



Having calculated the energy imparted to the screen, the next aspect of the invention 
concerns the removal from this energy value of the contribution to it by intensifier screen 
glare m step (4) which is directly related to hght exposure to the fihn. In summary this is 
done using a point spread function calculated by assuming that the x-ray photons are 
absorbed equally across each part of the intensifying screen corresponding to one pixel and at 
different depths within the screen. Knowing the pixel size, the solid angles between each 
potential absorption site and the neighbouring pixels can be computed. These solid angles 
define the proportion of light photons that the neighbouring pixels receive. The solid angles 
are weighted by distance from the absorption site to the pixel (representing light absorption 
by the screen) and by the actual x-ray energy reaching that site. Using these values glare can 
be estimated and removed from an image. The mathematics of this will now be explained in 
detail. 
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To derive the point-spread fiinction, the screen is considered to be divided into layers 
and then each layer to be sub-divided into sub-pixel size units. Each of these units is 
considered to be a potential site of absorption of an x-ray. Let t, cm be the thickness of the 
intensifying screen. The screen is considered to be in n layers so that each layer has a 
thickness of dt= t/n cm. The layer is subdivided into pixels (whatever resolution is being 
used) and then each pixel is split into 100 smaller elements. It is assumed that the x-ray film 
lies directly on top of the intensifying screen separated only by a screen overcoat so that there 
is a gap of /<,cm between the x-ray film and the intensifying screen. For each layer a 
weighting mask w/x, y) is computed (effectively the point spread function for that layer) 
which gives the percentage of light photons emitted, at (x„ yj, depth z, that reach the film 
corresponding to the spatial position (x, y). The proportion of photons from (x^ yJ reaching 
(x, y) is related to the solid angle dG from (x^ yJ, depth z to Cx. y). Symmetry round the 
azimuthal angle is assumed and the ID case considered (see Figure 32): 



15 de= CO A (•^2-^,)'-^'-c^' 



20 



Where b. c are found from simple geometry and x, and x, are the distances to either 
side of the target pixel (x, y) from the central pixel (x„ yJ. 

Some of the light photons that are emitted from (x^ y^ z) are absorbed by the 
phosphor so that de is weighted using Beer's law and the relative glare becomes: 



relative glare = dde 

light 



^ere Mplosphor is an average linear attenuation value. Dividing the relative glare by 
the total relative glare gives a weighting mask for each layer. The x-ray energy being 
25 imparted to each layer is now incorporated. The energy into each layer is: 
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where E.-. is the energy in,par,ed and is *e mcident energy. So .he relative giare now 

becomes: 



^phosphor ~ 



relative glare = dOe El"'''{x ,y ) 

These values are combined for each layer and for each sub-pixel to get the full 
5 weighting mask y) and the results are scaled so that yj=, .o and rotated around 

the central pixel to give a two dimensional mask. 

The parameters z/^*^ • i j 

^phosphor . Mphosphor ^rc mcludcd in the model more for illustrative 

purposes than practical value since it turns out that they have little effect on the PSF which 
IS the cuicial factor. Figures 26 and 27 show the modulation transfer functions (MTFs) of 
10 screens in which the varying parameters are and The modulation 

transfer function is simply the Fourier transform of the point spread function. Figure 26 
shows the MTF of intensifying screens using the model with 100 micron screen thickness 
and 0 gap between the film and the screen. There are three curves representing three 
different values of ^^Uor • As the value rises the more direct path to the film 
1 5 becomes more heavily favoured and so the MIT improves. Figure 27 shows the MTF of 

intensifying screens using the model with 100 micron screen thickness and 0 gap between the 
film and the screen. There are three curves representing three different values of 

''^phosphor ' 

As the value rises more x-r^ys get absorbed nearer the fihn and so the MTF impmves 

The two most important parameters in the model are the screen thickness , and the 
20 gap between the Hhn and the screen r. Figure 28 shows the PSFs for nvo differem screen 
thtclcnesses. The graphs show d>e percentage of light photons emitted from the central pixel 
(where the peak is ) and travelling to the surrounding pixels. The graphs are p^files across 
the 2D masks for when the spatial resolution is 50 microns. On the left, the screen is 100 
microns thick, whilst on the right its 200 microns thick. The thicker the screen a,e greater the 
abson„ion of the x-ray photons at depths furter away from the film and thus the more spread 
out the light emitted by the screen. The MTFs (for zen> gap) are shown in Figure 29 Figures 
30 and 31 show the PSFs and MTFs of when the gap between the fihn and the screen is 
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varied. Figure 30 shows the PSF of two intensifying screens and using the model with 100 
micron screen thickness and 0 and 100 micron gap between the film and the screen as 
marked. Figure 3 1 shows the MTF of two intensifying screens using the model with 100 
micron screen thickness and 0 and 100 micron gap between the film and the screen as 
marked. Clearly, these parameters have a large effect on the screen performance and they are 
an important part of the calibration data. For example the intensifying screen Fiji UM 
Mammo Fine enclosed within a Fuji EC-MA cassette has a thickness of approximately 
100 microns. The gap between film and screen is much more difficult to estimate and 
infonnation from manufacturers is sparse. However, there is probably at least a 1 micron 
protective coating on the film emulsion and a further protective coating of 5-1 5 microns for 
the screen. 

To properly compensate for glare it is necessary to know exactly where the edges of 
the mammographic film are since no glare comes from outside the film and those areas 
should be treated as zero in the convolution. Also, the energies imparted outside the breast 
area but on the film saturate the fihn so that to model the true effect of glare the energies on 
those regions have to be set to the expected incident energy as computed using the known 
time-of-exposure. The weighting mask w defines the point spread function for the 
intensifying screen so that the energy imparted image attained is simply the result of the 
energy imparted without glare convolved with w: 

20 ^-glare _ ^-no glare 

This can be solved by deconvolution in the Fourier domain. Figure 8 shows an 
example of the effects of compensating for the glare. 



15 



25 



30 



Noise detection 

Although the glare has been removed to enhance the calculation of /i,„, in steps (5) to 
(9) above, this step is itself important in that it allows the detection of microcalcifications and 
in particular allows a distinction to be made between microcalcifications and film-screen 
"shot" noise which can look confusingly similar to an image analysis system. 

Fihn-screen "shot" noise can arise from dust and dirt on the intensifying screen or 
from deficiencies in the film. A major difficulty in detecting microcalcifications is that this 
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no,se .ends ,o appear „i* si„,i,ar cha,ac«ris.ics to fte calcification: small, low filn, density 
bnght) and high frequency. Consequently, automated detection of micn,calcincaUons as 
locahsed bright spots tends to generate many faise positives. Although some of these can be 
elm,ma.ed by using the clustering property of real microcalcificaUons it would be preferable 
«> be able to eliminate them indiv.dually. The present invention uses the absence of blur from 
bnght spots to mark them as noise. The absence of blur indicates they were introduced into 
the nuagtng cham after the glare ftom the intensifying screen. 

In the above method of glare removal the mammogram is fust transfonned into an 
.mage which represents the energy imparted to the intensifying screen, h, this representation 
ca cficattons and noise have vety low values due to their apparent "high x-ray attenuation" 
When glare is removed from the image, pixels with low energy values due to noise be^me 
negative indicating that «>e original energy value was not feasible. Figure 9 shows two fito 
density profiles with confumed shot noise in them. It shows two profiles ac„ss ,00 micron 
.mages containing microcalcifications. The pixel value on the y-axis is linearly related to film 
dens.ty. with high pixel value meaning low fihn density and thus low energy imparted to the 
.ntensifying screen, h, the left profile the image restoration algorithm for glare picks out :^50 
as betng film-screen noise, and in the right profile the point at ;^206. The peaks at x-l 90 
and :c250 in the right profile are real calcifications. Figure 10 shows an original energy 
profile and the energies after glare compensation It shows two profiles ac„,ss a ,00 micton 
energy image. The left profile shows the enet^ies imparied before glare compensation TTtere 
■s a low value at ^2 which is due to noise (low energy means it appea. white). -Hte right 
pn^file shows the.energies after glare compensation, the energy at the noise point is now 
negative, i.e. i, has been mtroduced to the image after the major blurring stage 

Ftgure 1 1 shows an example of a physical phantom (TOR-MAX, University of Leeds 
Department of Medici Physics. X-ray Test Object) and the noise that the detection scheme ' 
finds. The image on the left is of a physical mammog„phic phantom. The phantom has 
stmulated microcalcifications (pointed at by the a.ow). round cylinders of varying densities 
and, at the far right, patches of texture. Before the mammogram was perfonned the 
film-screen combination was opened up and dust and dirt distributed randomly - these 
show-up as being bright white on the left image. The results of the noise-detection algorithm 
are shown on the right. All dte obvious noise pixels have been detected as well as some far 
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more subtle ones, none of the simulated calcifications have been marked as noise. 

Extensive testing of the film-screen shot noise detection algorithm has been cairied 
out on real mammograms. On a set of 20 sections of mammograms many of which contained 
microcalcifications an experienced radiologist marked the points corresponding to 
film-screen shot noise. She marked 156 points of which the above algorithm detected 150. 
The algorithm also detected 6 more points which the radiologist was unable to state 
categorically whether they were noise or real microcalcifications. No definite 
microcalcifications were marked as noise. Figuresl2 and 13 show examples of mammograms 
and the noise detected. Referring to Figure 12 on the left is a section of a real mammogram 
digitised to 1 00 microns per pixel. There is a cluster of noise points as confirmed by an 
experienced radiologist and by the algorithm as seen on the right which shows in white the 
noise detected. There is a quite obvious hair in the top right comer and some subtle genuine 
microcalcifications in the bottom left comer which have not been marked as noise. Referring 
to Figure 13 on the left is a section of a real mammogram digitised to 50 microns per pixel. 
There is a whole cluster of noise points as confimied by an experienced radiologists and by 
the algorithm as seen on the right which shows in white the noise detected. 

Although this noise detection scheme works extremely well, there are two 
considerations which should be noted: 

(i) At points where the film is saturated, that is where it has very high film 
density, or at very low film density such as in those areas beneath lead markers the estimated 
energy imparted is inaccurate and that can affect the noise detection. However, this happens 
only at the very edges of the breast image well away from any likely calcifications. An 
example of this is in Figure 1 1 where the very low film densities at the top of the image are 
marked as noise. 

points affect the noise detection at points near to them. The glare 
removal requires the energies in a local neighbourhood. If one of the energies is artificially 
low due to noise then the glare removal might be incorrect and might cause other noise not to 
be detected. This can be seen in Figurel 1 where there is a hair which clearly has the end 
points marked but not the central ones. This can necessitate a two-pass glare removal 
algorithm or require that any pixels nominated as calcifications within a certain radius of a 
noise point are also considered noise. 
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This med,od for de,ecd„g microcalcif.ca.i.ns can be furt,=r improved. „ „,ay be 

nite:/""™^""^^""'""— ^ 

.ha, U,e breas, co^.s. enUre.y of fa. and "i.,ere..i„g sh,ce calcifications have an 

-----o„coemcie„..h.ch.abo....™esh,.he..anU.os=«ss„e,,pes ^ 

an_ofanx-.,bea™«,™,ha„..ca,.fica.ion,pe.h^ 

comparable ,oU.a.U,rougb .3 .n, of i„«.es.ingdssue. For .his reason fte* value ' 

— aforpi.els„hichinfac.co.respond.oc..^^^^^^^^^ 

Which correspond to non-calcifications. 

Based on *is U,e improved n,eU,od uses ,he ..interesting tissue volume... denoted by 

Tcln 'bt " 

V,., can be computed f.„ the values over the region of interest on the image. Now 

cons.der the v,„ vaiue of a small volume of breast tissue, .hose achral volume is v., f no 
ca c,fica„on.s present in ir.v. Should be bound above by v... However, if ^ is a 
calcfication, ti,e computed v,„ would exceed v„ owing to the violation of tite fi,ndamental 
assump^onofthe.,., mode,. Thresholding the to v.„atio. enables differentiatioTof 
calctficattons from other breast tissue. This tatio is subject neither to varying imaging 
con ...ons nor to different tissue backgrounds on which a calcification is projected This is 
contia., ,0 .mage conttas, which ,s what most oti,er calcification detection algorithms 

con. '"""^ "'"^ --fy a weak 

congas constiam. and which are not too large. The calcification, of course, has oti.e bre^f 

ussueaboveandbelowitwhichbecauseoftheprojectivenantreofthemammogram 
~t. to the image. This contribution is removed by estimating the bac Jund 

cal ula ted. The estimation of v„ is more difficult, and to do ti,is the heuristic assumption is 
n,ad= that the candidate microcalcification has an elliptical cross section, so that its 
can be estimated from its projection in tiie image. 

In an initial experimental study a total of 20 image samples taken ftom 7 differem 
mammograms were used. The images are digitised to a resolution of 50m„ per pixel T^ere 
are altogether 27 microcalcifications in the 20 samples, each of which has at leas. , 

..age When ,hev.„,».hresho,dis se.,o3. The.go..hm detects.other re ions whiLe 
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radiologists are unable to make conclusive remarks on whether they correspond to real 
calcifications or non-calcifications. 

The results were also compared with those obtained by simply thresholding the grey 
level contrast. An ROC analysis shows that this method achieves both higher sensitivity and 
5 better specificity, albeit on a small sample. Figure 14 shows an example which compares a 
typical contrast-based detector attempting to detect a low-contrast microcalcification with the 
improved method above; 

(a) the original image with 2 microcalcifications; 

(b) detection result using the volume ratio algorithm: both calcifications are 
10 detected successftiUy; 

(c) detection result using grey level contrast: the more obscure calcification on the 
right is missed. 

Thus the glare removal process of the invention enables the elimination of many false 
positives that correspond to noise. Figure 15 shows an image sample containing image shot 
15 noise: 

(a) the original image with 1 microcalcification and a number of image shot noise 
points scattered over the image, of which the one lying horizontally near the bottom of the 
image is very prominent; 

(b) detection result using the volume ratio algorithm: the region coiresponding 
20 to shot noise is eliminated and only the real calcification is marked; 

(c) detection result using grey level contrast: the image shot noise is wrongly 
marked as a calcification. 

Algorithms based on contrast measures have difficulty rejecting such false positives. 
Thus despite the h,„, representation explicitly calcifications it is useful for detecting them and 
25 for differentiating false-positives due to shot noise. 



Anode heel eflect 



Returning to improvements in the calculation ofh,„. , the next step is to correct for the 
30 anode heel effect in step (5). 

The primary component of incident energy at (x. y) is directly proportional to the 
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number of photons incident to the volume of tissue projected onto that pixel: 

ETi^^y) = m.x,y)t^A^E';^Jx,y) , (3) 

where the last term is used to denote that part of Equation 2 which is independent of the total 
number of photons (nd stands for not-dependent). It is assumed that the x-ray energy 
5 spectrum stays the same but that the total number of photons ^ (V. x. changes with y) 
due to the anode heel effect and diverging beam. The scatter component at the pixel (. y) 
mostly comes from the x-ray photons that are entering the breast tissue in the surrounding 
neighbourhood. This neighbourhood is small enough to allow the anode heel effect over it to 
be Ignored so that the scatter component is also directly proportional to the incident radiation 
10 a.t(x,y): 

Erix,y) = m,x,yXA^Er„,{x,y) , (4) 

The total energy imparted is the sum of the primary and scatter components so that 

using Equations 3 and 4 gives: 

^'"^"'^^ = ^^^'^'^X^.(^r„/->;^) + ^r„,(:c,3.)). (5) 
1 5 The incident photon flux is greatest underneath the anode, let the position on the fihn 

at this pomt be (.^ yj. The aim is to change E'^''(x, y) to be as if from that incident photon 
flux. 

correcreAX,y) - .(^^^) ^ 

The ratio of the two photon fluxes needs to be computed. To detennine the anode 
20 heel effect for a specific system and to thus compute the ratio an x-ray exposure is performed 
wxth no object present which gives an apparently "blank fihn". For example, nine points on 
the film can be sampled. The energy imparted to the screen comes mostly from the primary 
radiation smce there is no scattering material: 

EZ.k(.x,y) = KV„x,yy^A^{E^''^^(^x,y) + Er„^(x,y)) 
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Noting that since the parts of the signal not dependent on the number of photons are 
equivalent (there is no object) gives: 

Substituting this into Equation 6 allows compensation for the anode heel effect. 
Estimation of compressed breast thickness 

An important step in this method is the estimation of the compressed breast thickness 
H. This is needed for an accurate calculation of but is also useful in checking the 
accumulated radiation dose per unit volume, and in locating the nipple which is important for 
diagnosis. 

Many recent mammography systems have built-in analogue or digital thickness 
meters but their accuracy and precision are currently wanting. In clinical practice, most 
existing systems do not have such indicators. 

This aspect of the invention provides a robust and accurate method for estimating the 
compressed breast thickness from a mammogram using image processing and modelling 
techniques above. The estimation is based upon the existence of the "breast edge", a 
fatty area around each breast where the breast thickness steadily reduces to zero. Determining 
that area using image processing provides enough data to estimate the breast thickness when 
calibration data such as the tube voltage and exposure time are known. The technique can 
also be applied to mammograms taken previously. 

During mammography, the breast is compressed between two supposedly parallel flat 
compression plates. The compression causes the breast to spread out, so that over most of the 
plate the breast is of equal thickness. However, towards the edge of the breast the breast 
bulges like a balloon and there is a not a straight vertical edge; Figure 16 shows this 
schematically. This bulge is what we term the breast edge. 

The breast tissue is enveloped in two layers of fibrous tissue, the deep layer overlying 
the muscle, and the very thin superficial layer beneath the skin. The superficial layer is 
separated from the skin by 0.5 to 2.5cm of subcutaneous fat or areolar tissue. Joining the 
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aye. ,o .he are f.„e fibrous (Coopers .iga„,e„., Fro. U,is, i. is reasonable 

o assuo,e «,a. wha. is ,er„ed bere U,e breas. edge co.is. e„.,re,y of fa, .bough in ,oca,i.ed 
reg,ons, pa„,cular,y near fte nipple and near U,e ligaments, m will no. be s.ricUy ttue 
F,gure , 7 ,s a schematic of a cranio-cauda, mammogram (no. .o scale)^ The left pi«ure 
shows .he breas., while d,e righ, picre is .he mammogram. The shaded areas are mean, .o 
represent volumes of i„.eres.ing .issue. The arc deno.ed by .he le„er E lies righ. on U,e edge 
of .he breas. as seen from .he x-.y source. ti,a. is. an x-ray from .he source .angen.ially 
.ouches *e breas. before reaching E. There is very linle auenuarion of *e x-ray beam 
anywhere along E. Along .he arc de„o.ed by .he lener D. .here is a grearer .hicl^ess of .issue 
0 han a. E. bu. stil, rela.ively linle a..„ua.ion, since this is .he breas. edge and mos. of .he 
reas. tissue is fa.. Along ti.e arc C, «,ere is H cm of fa. so ti,a. ti,e x-ray a„enua.ion is ,ui.e 
arge bu. .. .s still unifonn along ti,e arc. This changes for arcs A and B. where ti,ere is a 
he.erogeneous mixtirre of .issues. The x-ray a«enua.ion along .hese arcs can va.^ from being 
due.o//cmof .o almost cm of in.eres.ing tissue. 

This method of detennining breas. thickness H is based-upon delimiting .he projeced 
breas. edge from dre imeriorof d,e breas, i.e .o de.en„ine .he arcs C and E ,n Figure .7 
The breast can be segmemed from d,e background .o de.ennine .he smooti, arc E on ti,e basis 
Of ti,e film dens...es. Then all ti>e pixels which are slighUy brighrer ti,a„ «,ose on E are 
marked .o ge, ano.her smooti. curve, since .his is s.ill in the homogeneous b.«.. edge This 
may be arc D. Eventirally. arc C is reached, which is still smooth, since i. is jus, witiun the 
homogeneous breas. edge region. ,f pixels are continued to be marked which are slightiy less 
dark. .. starts to mark pixels witi.i„ .he imerior of tite breas. These could be anywhere 
and numerous since .he in«Hor is a heterogeneous mix of tissues. Thus ti,e arc d,a. delunits 
^.n enor o ti,e breast (a. C) is tite exti-emum of the smoo* lso-intensi.y curves. From ti.is 

mode, . .s also c e^that the breasttissue between arccandti^ex-raysourcecomprises;^ 
cm of fa. Thus, .f C can be found, and given agood model of tite mammogmphic imaging 
process, the p,xel values along C can be used to estimate tiae amount of fa. which must have 
been be^een the pixels along C and the x-ray source thus giving the hreas. thickness H 

^'■^--='«°P°-ble approaches to finding the curve which delimits fte breast 

interior and then estimating the breast thickne<!<: n^. * tr ^ ^ 

^ ^^"'•^""^"^"e's to find the smooth curve by using the 

ongmal pixels values which are produced by the dieitiser Th. • • 

cu oy me aigitiser. The imaging process is then 
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simulated using different breast thicknesses until the predicted pixel value matches that which 
.s found on the smooth curve. The second approach is to take a rough initial estimate of the 
breast thickness and then to generate the representation. Analysis of the h,„, representation 
enables determination of how the estimate of the breast thickness should be changed. When 
5 the A,, values are generated using an accurate breast thickness, it is found that within the 
projected breast edge there is so little attenuation that there cannot even be H cm of fat within 
that region. For these regions is set to zero and then it is determined what thickness of fat 
alone would give the observed attenuation. FigurelS shows some examples of the breast 
edges found by marking those pixels that have /i,, =0. The bright white areas represent the 
10 so-called breast edge. This is where the breast starts to decrease from its constant thickness H 
and the attenuation is no longer enough even to be ofH cm of fat. The inner edge of the 
projected breast edge region is mathematically quite smooth. Compensation for compression 
plate slant (typically up to 0.5cm from the chest-wall out to the nipple) is perfonned by 
effectively adding a wedge of fat to the breast by adjusting the mammographic image using 
1 5 monoenergetic modelling assumptions. 

The method of determining H is to start with an underestimate of the breast thickness 
and then to compute the A,,, representation. Those pixels that have h,^=H. h,=0 are marked 
and a measure of how rough the curve that those pixels represent is computed (see below). 
Initially, the h,„, values will be too high and there will be no pixels with h,„=0. As H 
20 increases, more and more pixels are detected with h,„,=0 and these pixels will represent a 
smooth curve. This continues until the internal breast region is reached at which point the 
roughness measure rises dramatically to indicate a rougher curve. 

Performing the computation this way rather than just using the pixel values in the 
original image and looking for iso-intensity curves allows more accuracy since H can be used 
25 to predict scattered and extra-focal radiation. It also enables checking that the algorithm is 
working properly. For example starting with many points that have > H then the value 
of// is far too low and should be increased it dramatically rather than in small steps. 
Another example is that if more that 20o/o of the breast is found as being within the projected 
breast edge, then //has increased too far. Figure 19 shows real examples of the breast 
30 edge as //varies. This shows images with luminance proportional to except for the bright 
white breast edge which is where A,, =0. The breast thicknesses tried are marked above each 
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■mage. The fi., value. 3.4c„, is fa. .00 ,ow - .here is no breas, edge; .he second value 5 4c™ 
.s JUS. abou. righ, - U,e.e is a suhs.an«ai breas, edge and a smooth in.ema, edge; ,he '.wo 
esumates are far ,00 high - Ae projected breas. edge has becon,e ragged 

As for the co.pu.a.io„ of roughness, as breas. thickness increases a largejump in *e 
measure rs expeced Jus. after U.e corvee, value. In Figure 20. *e curves marked A-E would 
Ideally have the roughness measures as shown on the right 

There are several possibiliries for such a roughness measure bu. U,e approach in the 
.nven.,on ,s ,0 use an es.imate of a fractal char^Cerisdc of 4e curve. A frac, curve has two 
paramere,.: n.c.a, dimension and D-dimension. and for es.ima«ng image .e«ures Ure ,a„er 
. more reliable, s^ble and gives better discrimination. A technique to estimate the 
D-dunension: "the covering blanket" which is based upon mon,ho,ogica. operaUons is used 
To compute the measure, two ftrther curves are creared from the curve: one from 
openmg the /,,„-« cun,e and one from closing it (see Fi^re 2,). An opening operation 
creates a smaller curve which is smooAer «,an *e original, whils. a closing operadon crea.es 
alargercut^ewhichissmootherthantheoriginal. ll.e area contained benveen *e .wo curves 
■s a measure of roughness: .he larger U,e area the rougher the curve. A 5 pixel by 5 pi.el disk 
■s used as our mon,hologica, operarcr and .he .hickness es.imadon is performed on 300 
micron resolution images. 

AIU,ough the area be^veen the curves is related to roughness, it is also direcUy rela.ed 
.0 the total leng* of Ote curve, so it is normalised and made dimensionless by dividing by 
the length of U,e curve. Taking all U,ese considera.ions into account the toughness is the 
number of pixels enclosed by closing and opening of divided by the number of pixels 
on the,=o curve. When .bis measure increases above a Areshold of 1.3 increase of/, is 
Stopped. 

In this es,ima.e of breas. Aickness. the smood, curve which delimirs .he imerior 
breas. region is de.ec.ed. The presence of Coopert ligamenrs in .he image, which is ,ui.e 
unusual, would afTec. only a secion of the curve and only slightly the over.,, roughness 
measure. The increased density of tissue around the nipple Jus. means that the curve is no. a 
sem,-c,rcle bu. a more irregular shape .ha. is s.i„ smoo*. Indeed, the fact Uta. .he curve is 
less regular can be used « de.ec. the posi.ion of Ae nipple 

The,, are several ways of fonning an initial (under)estima,e of breast dtickness from 
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the calibration data and image. One way exploits the fact that near the chest wall there is low 
scatter and the breast tends to be fatty so one can assume h,„,=0. //^«= H. From the energy 
imparted H can be estimated. Another way is to assume that near the breast edge there is 
pure fat and some nominal, high scatter-to-primary ratio. Another way is to estimate an 
5 initial value for the breast thickness using a film density that is known to be outside the breast 
and the calibration data; this always gives an underestimate of the actual breast thickness. 

Bounds on the breast thickness can also be estimated to check that the method is not 
trying infeasible breast thicknesses. The lower bound on H is related to the minimum 
attenuation apparent within the breast image. To achieve such low attenuation requires a 
10 certain minimum thickness of breast tissue. The minimum possible H occurs if the breast 
tissue has very low fat content so that /;,„,= H (i.e. highly attenuating). An upper bound on H 
can be determined in exactly the same way except using the maximum attenuation and 
considering the breast to be nearly all fat. 

As indicated in the introduction, the scattered radiation can be estimated using a 
1 5 published technique, achieving step (7), and then the extra-focal radiation must be removed in 
step (8). 

In this embodiment the extra-focal radiation is found by using the "inner projected 
breast edge" found in the estimation of// above. 

The extra-focal radiation is assumed to be travelling in random directions much the 
20 same as light from the sun on a cloudy day and is arriving at a point on the intensifying 

screen from all directions equally. If the thickness and the composition of the breast along 
some curve C is known, the expected primary and scatter components along C can be 
estimated. The extra-focal component along C can be found by subtracting them from the 
known total energy imparted; 

25 ET' = E';^ - ET" - 

this value is then extrapolated over the entire image. 

As mentioned above the curve C used is the "inner projected breast edge". This is the 
breast tissue around the edge of the compressed breast which is fat and which is where the 
breast starts to curve. Figure 3 shows this pictorially. In Figure 3 it can be seen that most 
30 x-ray radiation comes from the focal spot. However, some extra-focal radiation is scattered so 
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15 



.ha, „ reaches U,e fihn-screen and co„,ribu.=s ,o U,e fmal image. Since the breas. edoe is 
cun,ed, photons IB arriving at low angles have easy access to the film-screen and thus have 
a stgntficant impact on the final image. In the following the breast edge is defined to be that 
volume of the breast which is fat and is where the breast thickness descends fh,m H towards 
zero. On a mammogram this area can be depicted by an inner and outer curve 3 A and 3B as 
shown m d,e top-righ, of flgure 3 and these can be calculated using d,e techniques described 
above for the calculation of H. 

The inner curve is used as curve C and «,e breast edge is modelled as consisting 
entirely of fat and being semi-circular. Using this assumpUon. the thicbtess of tissue 

1 0 between any point 

and where the e>c,ra-focal is perceived to be coming from can be computed and then this is 
used to adjust the extta-focal estimate accordingly. 

Figure 22 shows the coordinate system used. (X. Y. Z) are used as the 3D coordinates 
of a pom, with the JTaxis along the chest wall, ^ Z axis perpendicular to the chest wall and 
Y the vertical axis towards the x-ray source, ft is a point in the 2D image plane. When no 
breast is present the film-screen combination at ft z) is assumed u, receive the same 
extra-focal radiation from all directions ( 6. p). p is termed the azimuthal angle 6 the 
elevation angle and is .be "extm-focal constant", that is, the energy incidem on ft 

from any direction. 

The energy imparted to the screen due to extra-focal radiation can be assumed to be 
constant across the image if there were no breast present. With the breast present, the 
extra-focal radiation is attenuated by that quantity of breast tissue present along any particular 

path of travel. 

That quantity of breast tissue is estimated using a model of the breast edge and the 
25 constant that would have been imparted had no breast been present is computed 

At each pixel (x, z). the energy imparted when a grid is used can be approximated by: 

'0 where T is the transmission through the anti-scatter grid, and assuming that the screen 
absorbs all the photons reaching it. regardless of angle. To simplify the analysis a 



20 
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monoenergetic case with photon energy e is treated. The total extra-focal radiation imparted 
to the intensifying screen at any point when no object is present is given by: 

2/T;r 2,-r 
0 0 0 ° 

5 Let /J3(x, z, 0, p) be the attenuation due to the breast along the path to point (x, z) from 

angles : 6, p. 

pKx,z,9,p) = K.{x,z,d,p)p.^^ie)^h^^{x,z,d,p)p^^X^) 

with s the relevant photon energy. The energy imparted due to extra-focal radiation can be 
estimated using Beei^s Law, ignoring scatter, for when the breast is present: 



10 Er{x,z) = E'z!\\le-^^--'-^^Tie,p)dedp 



To study the problem analytically, this equation must be simplified. This is done by 
first assuming symmetry around the azimuthal angle p so that the problem becomes 
essentially one dimensional: 

Er(x,z) = iTTElZ^j e-^^'-''^T(e)de 

0 

1 5 E™P„,„ is estimated together with the attenuation along each ray using a suitable 

model of the breast and breast shape as explained below. The problem can be reduced further 
since 7c/2<e <k represents the angles of extra-focal radiation coming from the chest wall 
where there is a full thickness of breast and body, so that contribution is assumed to be 
minimal. Furthermore, the most important part of the breast for extra-focal radiation is the 

20 breast edge and that can be considered to be just fat, so that the equation can be further 
simplified: 



2 



Er(x,z) = 2;rEZ4 e''^''^'''''''T{e^de , (7) 




10 



-30- 



In this equation h,Jx, z. 0), i.e. the thickness of fat from any point along any direction 
is computed as explained below. Using the coordinate system of Figure 22 the shape of the ' 
breast edge is considered to be semi-circular. With no loss of generality save for assumptions 
of rotational symmetry about the azimuth, the ID case where X=0 and Z is perpendicular to 
the chest-wall is treated, so that the breast edge can be described by: 

Y{Z,) = a±^ia'-Zl), 

where a =H/2 andZ^ is the Z-coordinate relative to the edge of the inner breast edge, i.e. Z, = 
Z - Z,,,,. Figure 23 shows the two cases inside the breast. The case Al shows the ray 
travelling through both breast and breast edge: for angles of 0 < 0 <e„ where 0„ = tan- H/| Z„ 
I the extra-focal rays must travel through a thickness of breast tissue and a thickness of breast 

edge tissue so that: 



cos e cos e 

and Z,„,_, is the intersection of the ray and the breast edge which occurs at: 

intersect ~ 

\ + nr 

1 5 where m = ^ and c = Z, tan 6. For the points which satisfy case A2, some of the 

extra-focal rays will come through tissue which is within the uniform breast thickness and 
some through the breast edge area. These rays are from e„ to ./2 and the thickness of the 
breast is simply: h.je) = Wsin 6. Outside the /,,,,=0 line, there are angles that have a free 
line-of-sight to the pixel (x, z), see case Bl in Figure 23. In fact, these angles are from 0 to 
eu„, where is the angle of the tangent to the breast curve which passes through the 
pixel. The tangent to the bottom half of the breast which passes through the point (Z, 0) that 
is required passes through the point: 



20 



and thus: 
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^umit = tan-' 
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z 

^tangent 



So, that for 0 <9<eii„^ : hfJO) = 0. Case B2 in Figure 23 reveals that for some of the 
angles of interest the x-ray path is purely through the breast edge. Now the two intersection 
points of the ray with the breast edge semi-circle are needed, i.e.(Zi„,„ect i, ^-^^^x i) and 

^ intellect 2> "^intersect 2)- 

The two Z coordinates are given by: 

m{a - c) ± ^(m^a^ - -h lac) 

with m -tan 0, a =H/2 and c = -tan OZq. And thus: 



^fat ) V^^'^intersect 1 ^intersect 2) (intersect 1 intersect 2) 

Using these thicknesses the basis extra-focal equation (7) can be rewritten as: 

W2 



\ T(^)d^+ jQ''^''^-^''''''^T(0)dd 

V 0 0 



(8) 



where 0^^^ is the angle of the tangent to the breast edge curve and for cases Al and A2 is 
15 zero. 

In order to finally compute the extra-focal radiation from equation (8) it is necessary 
to estimate the constant The location of the projected breast edge is assumed and 

since the breast is assumed to be all fat the location of a curve which has hi„r=0 and hyi,,=H 
20 (the inner curve in Figure 3). Thus the primary and scatter components along the curve can 
be estimated and by 
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subtracting those values away from the total energy the extra-focal component alone the 
curve can be found. Knowing just one value for the extra-focal component allows 
computation of the constant £'"'^„,,„ as follows: 

Let the curve with ;,,,,=0 and be C.,,=,. Since the location of that curve is 
5 known from the breast edge an estimate of the extra-focal component there is: 

The scatter estimate comes from step (7) of the method, „ comes from the image 
and comes from knowing that A,,,=0 and V=H. The average of along the C,„,-, 
curve and Equation (8) with 9 = 0 to detennine the constant: 



exlra Jl2 — 

2k J e-''^''''^''^'---'^T{e)dd 



Thus the extra-focal component and the percentage of the total radiation that is 
extra-focal can be computed. Figure 24 shows two examples of the extra-focal componem 
computed in this way. They show that, as expected, the extra-focal radiation is high at the 
breast edge but falls rapidly inside the breast falling to zero if the breast is large enough. 

Incoiporating all the assumptions above in examples, the percentages computed 
indicate that extra-focal radiation is of the order of 6-IO0/0 of the total radiation. This is on the 
low sxde of the estimate made in other ways which state that extra-focal component can make 
up to 15% of the total radiation, but is thought to be reasonable. 



Calculation of /i^^ 



The representation can now be found by applying steps (9) and (10) to produce h 
surfaces such as that shown in Figure 5. In practice, steps (6) to (10) are used to calculate a "" 
first value of and then these steps are repeated with improved assumptions about the 
breast thickness until the best value of A,„, is found. 



• • • 



CLAIMS 



1 . A method of correcting a digital representation of an x-ray image for 
5 degradation caused by a digitiser, the digital representation being a measurement of the image 
density for each of a plurality of pixels of the image, the measurement being obtained by 
illuminating the image, measuring the attenuation of the light by the image and calculating 
from the attenuated light values the image density, the method comprising the steps of: 

calculating from the pixel values in the digital representation the value of the intensity 
10 of said attenuated light corresponding to each pixel; 

applying to the intensity values a modulation transfer function to correct for the 
degradation introduced by the digitiser; and 

converting the corrected intensity values back into values representative of the image 
intensity. 

15 

2. A method according to claim 1, wherein the modulation transfer function 
corrects for blur and the digitising aperture. 

3. A method according to claim 1 or 2, wherein the intensity 7) of the attenuated 
20 light is calculated as: 7"/=// x 70"^, where Ij is the intensity of the illumination of the image and 

D is the value of the image density in the digital representation, and the corrected intensity 
values T/ are converted back into values D' representative of the image 
intensity by: D* = log(I/Ti 

25 4. Apparatus for correcting a digital representation of an x-ray image for 

degradation caused by a digitiser, the digital representation being a measurement of the image 
density for each of a plurality of pixels of the image, the measurement being obtained by 
illuminating the image, measuring the attenuation of the light by the image and calculating 
from the attenuated light values the image density, the apparatus comprising: 

30 calculating means for calculating from the pixel values in the digital representation 

the value of the intensity of said attenuated light corresponding to each pixel; 
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correction means for applying to the intensity values a modulation transfer function to 
correct for the degradation introduced by the digitiser; and 

conversion means for converting the coirected intensity values back into values 
representative of the image intensity. 

5. Apparatus according to claim 4. wherein the conversion means corrects for 
blur and the digitising aperture. 

6. Apparatus according to claim 3 or 4, wherein the calculation means calculates 
intensity T, of the attenuated light as: T.=I, x JO- where I. is the intensity of the illumination 
of the image and D is the value of the image density in the digital representation, and the 
conversion means converts the corrected intensity values T,' back into values D' 
representative of the image intensity by: D' = log(I/T, '). 

7. A computer program storage medium readable by a computer system and 
encoding a computer program for controlling a computer to correct a digital representation of 
an x-ray image for degradation caused by a digitiser, the digital representation being a 
measurement of the image density for each of a plurality of pixels of the image, the 
measurement being obtained by illuminating the image, measuring the attenuation of the light 
by the image and calculating from the attenuated light values the image density, by a method 
comprising the steps of: 

calculating from the pixel values in the digital representation the value of the intensity 
of said attenuated light corresponding to each pixel; 

applying to the intensity values a modulation transfer function to correct for the 
degradation introduced by the digitiser; and 

converting the corrected intensity values back into values representative of the image 

intensity. 



8. A computer program storage medium according to claim 7, wherein the 
modulation transfer function corrects for blur and the digitising aperture. 
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9. A computer program storage medium according to claim 7 or 8, wherein the 
intensity 7/ of the attenuated light is calculated as: Ti^I^x IQr^^ where \ is the intensity of the 
illumination of the image and D is the value of the image density in the digital representation, 
and the corrected intensity values are converted back into values U representative of the 

5 image intensity by: Z) ' = log(I/Ti 

10. In a method of x-ray imaging using an intensifying screen to receive x-rays 
and emit light to be recorded on an x-ray film, a method of calculating from the fihn image 
density the energy which was imparted to the intensifying screen, comprising the steps of: 

1 0 calibrating the film and intensifying screen by measuring the response of the film and 

intensifying screen to a plurality of different intensities of received x-rays; 

fitting a theoretical model of the expected response to the measured response, and 
using the fitted theoretical model to calculate the imparted energy from the image density. 

1 ^ 11- A method according to claim 1 0, wherein the theoretical model is a serpentine 

curve of the form: 

20 where a and b are constants, x is the logarithm of the energy imparted to the 

intensifying screen and y is the image density. 

12. A method according to claim 10 or 1 1, wherein the response of the film and 
intensifying screen to a plurality of different intensities of received x-rays is measured by 

25 exposing the film-screen combination to x-rays through a lucite step wedge and measuring 
the image density produced by the exposure through the different steps of the wedge. 

13. Apparatus for calculating from image densities of an x-ray film image the 
energy which was imparted to an intensifying screen used to receive x-rays and emit light to 

30 be recorded on the x-ray film, the apparatus comprising:- 

fitting means for fitting a theoretical model of the expected response of the x-ray film 
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and intensifying screen u. a plurality of different intensities of received x-.ys to a measured 

response, and 

calculation means for using the fitted theoretical model to calculate the imparted 
energy fi-om the image density. 

5 

14. Apparatus according to claim 13, wherein the theoretical model is a serpentine 
curve of the form: 
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where a and 3 are constants, x is the logarithm of the energy imparted to the 
intensifying screen and y is the image density. 

15. A computer program storage medium readable by a computer system and 
encoding a computer program for controlling a computer to calculate from image densities of 
an x-ray film the energy which was imparted to an intensifying screen used to receive x-rays 
and emit light recorded on the x-ray film, by a method comprising the steps of 

fitting a theoretical model of the expected response of the x-ray film and intensifying 
screen to a plurality of different intensities of received x-rays to a measured response and 
using the fitted theoretical model to calculate the imparted energy from the image density. 

1 6. A computer program storage medium according to claim 1 5, wherein the 
theoretical model is a serpentine curve of the form: 

x^y+ a^y-b^x= 0, 

where a and ^, are constants, . is the logarithm of the energy imparted to the 
intensifying screen and y is the image density. 



17. In a method of x-ray imaging in which an intensifying screen is 



used to 
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receive x-rays and emit light to be recorded on an x-ray, the image recorded on the film is 
digitised to give a pixelised digital representation of the image density on the film, a method 
of enhancing the digital representation of the image to remove the contribution thereto of 
glare from the intensifying screen comprising the steps of 

converting the digital representation of image density into a representation of the 
energy imparted to the screen; and 

deconvolving the representation of the energy imparted to the screen using a 
weighting mask defining the point spread function for the intensifying screen to remove the 
contribution thereto of glare from the intensifying screen. 

18. A method according to claim 1 7, wherein the weighting mask defines the 
contribution to the energy imparted to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels weighted by distance from the 
absorption site to the pixel. 

1 9. A method according to claim 1 8 , wherein the weighting mask defmes the 
contribution to the energy imparted to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels further weighted by the depth of 
the absorption site in the intensifying screen. 

20. A method of detecting fihn-screen shot noise in a x-ray image in which an 
intensifying screen is used to receive x-rays and to emit light to be recorded as an image on a 
x-ray film, the image recorded on the film being digitised to give a pixelised digital 
representation thereof which is then converted into a representation of the energy imparted to 
the screen, the method comprising removmg from that representation glare from the 
intensifying screen in accordance with the method of claim 17. and defining as fihn-screen 
shot noise regions of the image in which the resultant value for energy imparted to the screen 
is infeasible. 

21. A method according to claim 20, wherein regions of the image in which 
resultant values for energy imparted to the screen having negative values are defined as film- 



screen shot noise. 
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22. A method according to claim 20 or 2 1 . wherein the 



image is a mammogram. 



23. Apparatus for enhancing the digital representation of an x-ray image obtained 
by using an intensifying screen to receive x-rays and emit light to be recorded on an x-ray 
film and digitising the film to give a pixelised digital representation of the image density 
thereon the apparatus comprising: 

conversion means for converting the digital representation of image density into a 
representation of the energy imparted to the screen; and 

deconvolution means for deconvolving the representation of the energy imparted to 
the screen using a weighting mask defining the point spread function for the intensifying 



screen to remove the contribution thereto of glare from the intensifying 



screen. 



24. Apparatus according to claim 23, wherein the weighting mask defines the 
contribution to the energy imparted to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels weighted by distance from the 

absorption site to the pixel. 



25. Apparatus according to claim 24 , wherein the weighting mask defines the 
contribution to the energy imparted to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels further weighted by the depth of 
the absorption site in the intensifying screen. 

26. Apparatus for detecting film-screen shot noise in a x-ray image obtained by 
using an intensifying screen to receive x-rays and to emit light to be recorded as an image on 
an x-ray fihn, the image recorded on the film being digitised to give a pixelised digital 
representation thereof which is then converted into a representation of the energy imparted to 
the screen, the apparatus comprising apparatus according to claim 23, 24 or 25 for removing 
from that representation glare from the intensifying screen, and means for defining as film- 
screen shot noise regions of the image in which the resultant value for energy imparted to the 
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screen is infeasible. 

27. Apparatus according to claim 26, wherein regions of the image in which 
resultant values for energy imparted to the screen having negative values are defined as film- 

5 screen shot noise. 

28. Apparatus according to claim 26 or 27, wherein the image is a mammogram. 

29. A computer program storage medium readable by a computer system and 

1 0 encoding a computer program for controlling a computer to enhance the digital representation 
of an x-ray image obtained by using an intensifying screen to receive x-rays and emit light to 
be recorded on an x-ray film, the image recorded on the film being digitised to give a 
pixelised digital representation of the image density thereon, the removal from the digital 
representation of the image of the contribution thereto of glare from the intensifying screen 

1 5 comprising the steps of: 

converting the digital representation of image density into a representation of the 
energy imparted to the screen; and 

deconvolving the representation of the energy imparted to the screen using a 
weighting mask defming the point spread function for the intensifying screen to remove the 
20 contribution thereto of glare from the intensifying screen. 

30. A computer program storage medium according to claim 29, wherein the 
weighting mask defines the contribution to the energy imparted to the screen caused by the 
emission of light from the point of absorption of an x-ray towards neighbouring pixels 

25 weighted by distance from the absorption site to the pixel. 

31. A computer program storage medium according to claim 30 , wherein the 
weighting mask defmes the contribution to the energy imparted to the screen caused by the 
emission of light from the point of absorption of an x-ray towards neighbouring pixels further 

30 weighted by the depth of the absorption site in the intensifying screen. 
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32. 



A computer program storage medium readable by a computer system and 
encodmg a computer program for controlling a computer to detect filn.-screen shot noise in a 
x-ray .mage obtained by using an intensifying screen to receive x-rays and to emit light to be 
recorded as an image on an x-ray film, the image recorded on the film being digitised to give 
a p.xehsed digital representation thereof which is then convened into a representation of the 
energy imparted to the screen, the detection of film-screen shot noise comprising removing 
from that representation glare from the intensifying screen in accordance with claim 29 Jd 
defining as film-screen shot noise regions of the image in which the resultant value for 
energy imparted to the screen is infeasible. 

33. A computer program storage medium according to claim 32, wherein regions 
of the .mage in which resultant values for energy imparted to the screen having negative 
values are defined as film-screen shot noise. 

15 34. A computer program storage medium according to claim 32 or 33. wherein the 

image is a mammogram. 



35. A method of calculating from a mammogram the compressed thickness of the 



20 



imaged breast, comprising the step of delimiting in the mammogram the region 
coxxesponding to the part of the breast which is compressed from the region corresponding to 
the uncompressed breast edge by detecting the smoothness of curves of equal intensity in the 



mammogram. 



25 
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.6. A method according to claim 35. further comprising the steps of detecting the 
smooth equal intensity cun^e which is positioned furthest from the breast edge, calculating 
from the mtensities along that curve the thickness of fat in the breast which gives those 
intensities, and calculating the thickness of the compressed breast from that thickness of fat. 

37. A method according to claim 35 or 36. comprising as an initial step the 
conversion of the mammogram into the /,,„, representation representing the thicknesses of fat 
and mteresting tissue in regions of the breast contributing to each point in the mammogram 
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38. A method according to claim 37, wherein the conversion is performed on the 
basis of an underestimate of the compressed breast thickness, the conversion being iterated 
with successively changing estimates of the compressed breast thickness until a smooth 
delimitation between the region corresponding to the part of the breast which is compressed 
and the region corresponding to the uncompressed breast edge is found. 

39. A method according to any one of claims 35 to 38, wherein a predefined 
smoothness threshold is set to detect the smooth curves. 

40. A method of calculating the contribution to a mammogram of extra-focal 
radiation by calculating the compressed thickness of the imaged breast according to the 
method of any one of claims 35 to 39, calculating from the x-ray exposure and an estimate of 
the scattering of radiation the expected intensity in the mammogram along the delimitation 
between the region corresponding to the part of the breast which is compressed and the region 
corresponding to the uncompressed breast edge, and comparing the expected intensity to the 
actual intensity in the mammogram. 

41. Apparatus for calculating from a mammogram the compressed thickness of 
the imaged breast, comprising means for delimiting in the mammogram the region 
corresponding to the part of the breast which is compressed from the region corresponding to 
the uncompressed breast edge by detecting the smoothness of curves of equal intensity in the 
mammogram. 

42. Apparatus according to claim 41, further comprising detection means to detect 
the smooth equal intensity curve which is positioned furthest from the breast edge, and 
calculating means for calculating from the intensities along that curve the thickness of fat in 
the breast which gives those intensities, and for calculating the thickness of the compressed 
breast from that thickness of fat. 



43. Apparatus according to claim 41 or 42, comprising conversion means for 
converting the mammogram into the representation representing the thicknesses of fat and 
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interesting tissue in regions of the breast contributing to each point in the mammogram. 

44. Apparatus according to claim 43, wherein the conversion means first 
underestimates the compressed breast thickness, and iterates the conversion with successively 

5 changing estunates of the compressed breast thickness until a smooth delimitation between 
the region corresponding to the part of the breast which is compressed and the region 
corresponding to the uncompressed breast edge is found. 

45. Apparatus according to any one of claims 41 to 44, wherein a predefined 
1 0 smoothness threshold is set to detect the smooth curves. 



46. Apparatus for calculating the contribution to a mammogram of extra-focal 
radiation comprising apparatus according to any one of claims 41 to 45 for calculating the 
compressed thickness of the imaged breast, and further comprising expected intensity 

1 5 calculating means for calculating from the x-ray exposure and an estimate of the scattering of 
radiation the expected intensity in the mammogram along the delimitation between the region 
corresponding to the part of the breast which is compressed and the region corresponding to 
the uncompressed breast edge, and comparing means for comparing the expected intensity to 
the actual intensity in the mammogram. 

20 

47. A computer program storage medium readable by a computer system and 
encoding a computer program for controlling a computer to calculate from a mammogram the 
compressed thickness of the imaged breast by a method comprising the step of delimiting in 
the mammogram the region corresponding to the part of the breast which is compressed from 

25 the region corresponding to the uncompressed breast edge by detecting the smoothness of 
curves of equal intensity in the mammogram. 

48. A computer program storage medium according to claim 47, wherein the 
method further comprises the steps of detecting the smooth equal intensity curve which is 

3 0 positioned furthest from the breast edge, calculating from the intensities along that curve the 
thickness of fat in the breast which gives those intensities, and calculating the thickness of the 
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compressed breast from that thickness of fat. 



49. A computer program storage medium according to claim 47 or 48, wherein the 
method further comprises as an initial step the conversion of the mammogram into the 
representation representing the thicknesses of fat and interesting tissue in regions of the breast 
contributing to each point in the mammogram. 

50. A computer program storage medium according to claim 49, wherein the 
conversion is performed on the basis of an underestimate of the compressed breast thickness, 
the conversion being iterated with successively changing estimates of the compressed breast 
thickness until a smooth delimitation between the region corresponding to the part of the 
breast which is compressed and the region corresponding to the uncompressed breast edge is 
found. 

51. A computer program storage medium according to any one of claims 47 to 50, 
wherein a predefined smoothness threshold is set to detect the smooth curves. 

52. A computer program storage medium readable by a computer system and 
encoding a computer program for controlling a computer to calculate the contribution to a 
mammogram of extra-focal radiation by calculating the compressed thickness of the imaged 
breast according to the method of any one of claims 35 to 40, calculating from the x-ray 
exposure and an estimate of the scattering of radiation the expected intensity in the 
mammogram along the delimitation between the region corresponding to the part of the 
breast which is compressed and the region corresponding to the uncompressed breast edge, 
and comparing the expected intensity to the actual intensity in the mammogram. 

53. A method of processing a mammogram to produce an representation 
representing the thicknesses of fat and interesting tissue in regions of the breast contributing 
to each point in the mammogram, the mammogram being produced by using an intensifying 
screen to receive x-rays transmitted through the breast and to emit light to be recorded as the 
mammogram on an x-ray film, the mammogram being digitised by a digitiser to obtain a 
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digital representation thereof comprising pixel values representing the image intensities in the 
mammogram, the method comprising the steps of; 

correcting the digital representation for degradation caused by the digitiser according 
to the method of any one of claims 1 to 3; 

converting the values in the digital representation to values representing the energy 
which was imparted to the intensifying screen according to the method of any one of claL 
10 to 12; 

compensating the converted values for glare from the intensifying screen according to 
the method of any one of claims 12 to 22; 

further compensating the values for variations in intensity of the x-ray beam incident 

on the breast; 

calculating the compressed thickness of the imaged breast according to the method of 
any one of claims 35 to 39; 

calculating the contribution to the mammogram of extra-focal radiation according to 
15 the method of claim 40; 

subtracting from the compensated values the extra-focal radiation and an estimate of 
the scattered radiation to give enhanced measured values of the energy which was imparted to 
the intensifying screen; and 

calculating the representation by equating the enhanced measured values with 
"0 predicted values. 
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54. A method according to claim 53, wherein the predicted values are calculated 
from the intensity of x-rays incident on the breast and the absoiption characteristics of 
physical elements of the imaging system. 

55. Apparatus for processing a mammogram to produce an h,„, representation 
representing the thicknesses of fat and interesting tissue in regions of the breast contributing 
to each point in the mammogram, the mammogram being produced by using an intensifying 
screen to receive x-rays transmitted through the breast and to emit light to be recorded as the 
mammogram on an x-ray fxlm, the mammogram being digitised by a digitiser to obtain a 
digital representation thereof comprising pixel values representing the image intensities in the 
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mammogram, the apparatus comprising: 

apparatus for correcting the digital representation for degradation caused by the 
digitiser according to any one of claims 4 to 6; 

apparatus for calculating from the values in the digital representation values 
5 representing the energy which was imparted to the intensifying screen according to claim 13 
or 14; 

apparatus for compensating the converted values for glare from the intensifying 
screen according to any one of claims 23 to 28; 

means for further compensating the values for variations in intensity of the x-ray 
10 beam incident on the breast; 

apparatus for calculating the compressed thickness of the imaged breast according to 
any one of claims 41 to 45 

apparatus for calculating the contribution to the mammogram of extra-focal radiation 
according to claim 46; 

1 5 subtraction means for subtracting from the compensated values the extra-focal 

radiation and an estimate of the scattered radiation to give enhanced measured values of the 
energy which was imparted to the intensifying screen; and 

hi^t calculating means for calculating the representation by equating the enhanced 
measured values with predicted values. 

20 

56. A computer program storage medium readable by a computer system and 
encoding a computer program for controlling a computer to process a mammogram to 
produce an hi^t representation representing the thicknesses of fat and interesting tissue in 
regions of the breast contributing to each point in the mammogram, the mammogram being 
25 produced by using an intensifying screen to receive x-rays transmitted through the breast and 
to emit light to be recorded as the mammogram on an x-ray film, the mammogram being 
digitised by a digitiser to obtain a digital representation thereof comprising pixel values 
representing the image intensities in the mammogram, by a method comprising the steps of: 

correcting the digital representation for degradation caused by the digitiser according 
30 to the method of any one of claims 1 to 3; 

converting the values in the digital representation to values representing the energy 
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which was imparted to the intensifying screen according to the method of any one of claims 
10 to 12; 

compensating the converted values for glare from the intensifying screen according to 
the method of any one of claims 12 to 22; 

further compensating the values for variations in intensity of the x-ray beam incident 
on the breast; 

calculating the compressed thickness of the imaged breast according to the method of 
any one of claims 35 to 39; 

calculating the contribution to the mammogram of extra-focal radiation according to 
10 the method of claim 40; 

subtracting from the compensated values the extra-focal radiation and an estimate of 
the scattered radiation to give enhanced measured values of the energy which was imparted to 
the intensifying screen; and 

calculating the h,„, representation by equating the enhanced measured values with 
15 predicted values. 



57. A method of detecting microcalcifications in a breast from a mammogram of 
the breast, comprising the steps of processing the mammogram to produce an h,„, 
representation according to the method of claim 53 or 54, converting the value ofh,„, for a 
candidate region into a value representative of the volume of interesting tissue and 
thresholding the volume values to detect as microcalcifications areas of the mammogram 
where the volume value exceeds a threshold. 
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58. A method according to claim 57 wherein the value of the volume of interesting 
25 tissue is compared to an estimate of the actual volume of the candidate region. 



59. A method according to claim 57 or 58 wherein the volume values 
calculated after subtraction from the A,„, value of an estimate of the background h. 



are 



30 60. Apparatus for detecting microcalcifications in a breast from a mammogram of 

the breast, comprising apparatus for processing the mammogram to produce an 
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ntation according to claim 55, means for converting the value of for a candidate 
nto a value representative of the volume of interesting tissue and means for 
Iding the volume values to detect as microcalcifications areas of the mammogram 
he volume value exceeds a threshold. 

6 1 . Apparatus according to claim 60 further comprising means for comparing the 
f the volume of interesting tissue to an estimate of the actual volume of the candidate 

62. Apparatus according to claim 60 or 61 further comprising means for 
ting from the value an estimate of the background /z/„„ before calculating the 
values. 

63. A computer program storage medium readable by a computer system and 
ig a computer program for controlling a computer to detect microcalcifications in a 
rom a mammogram of the breast by a method comprising the steps of processing the 
3gram to produce an /2/„^ representation according to the method of claim 53 or 54, 
ing the value of A/^, for a candidate region into a value representative of the volume of 
ing tissue and thresholding the volume values to detect as microcalcifications areas of 
nmogram where the volume value exceeds a threshold. 

64. A computer program storage medium according to claim 63 v^herein the value 
olume of interesting tissue is compared to an estimate of the actual volume of the 

ite region. 

65. A program storage medium according to claim 63 or 64 wherein the volume 
u-e calculated after subtraction from the value of an estimate of the background 

66. A method of processing a mammogram substantially as hereinbefore described 
erence to and as illustrated in the accompanying drawings. 
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Abstract 
X-ray image Processing 
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